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ly^ ' Abstract 

I In this paper we obtain a description of the Hermitian operators 

■ acting on the Hilbert space C", description which gives a complete 

solution to the over parameterization problem. More precisely we pro- 
vide an explicit parameterization of arbitrary rt-dimensional operators, 
operators that may be considered either as Hamiltonians, or density 
pH . matrices for finite-level quantum systems. It is shown that the spee- 

ch I tral multiplicities are encoded in a flag unitary matrix obtained as an 

^ ■ ordered product of special unitary matrices, each one generated by 

a complex n — /c-dimensional unit vector, fc = 0, 1, . . . , n — 2. As a 
byproduct, an alternative and simple parameterization of Stiefel and 
Grassmann manifolds is obtained. 



1 Introduction 



There is a considerable interest in a simple description of density matrices 
that have a wide variety of apphcations, particularly in quantum information 
theory, and many efforts were devoted in describing them. However the prob- 
lem of over parameterization is stih open JJj, for an arbitrary n. We solve 
this problem by providing an explicit parameterization of eigenvalues, as 
well as of the unitary matrices that diagonalize arbitrary finite-dimensional 
Hermitian operators. Quite expectedly, such a description is closely related 
to the description of various homogeneous manifolds, alluded to the title. On 
the other hand the Stiefel, or Grassmann manifolds arise in many problems 
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from different other domains such as encryption, coherent states, geometric 
phases, signal processing, geometric integration on homogeneous manifolds, 
numerical linear algebra algorithms, and many others. In such problems, 
the "states of interest" are, in general, elements of some homogeneous space 

X'^G/K (1.1) 

where G is a Lie group, and K isa. closed subgroup of G. Further in problems 
arising in engineering, physics, quantum information theory, one needs a 
concrete realization of these manifolds in a form able to be stored into a 
computer. Although, the geometrical description of Grassmann and Stiefel 
manifolds is available in many books, see for example 0, [Zj, the available 
parameterizations of the above manifolds are not the most convenient in 
some concrete applications, e.g. j3I, US] ■ 

One aim of the paper is to obtain a satisfactory description of Hermitian 
operators that appear in the study of finite-level systems. By satisfactory, we 
mean a complete as possible description of equivalence classes of Hermitian 
operators, i.e. a full and complete parameterization of the orbits gener- 
ated by these operators. We recall that two Hermitian operators are in the 
same orbit if their spectra coincide, or, equivalently, if their characteristic 
polynomials are identical. Our own interest in studying such problems was 
originally awakened by a query raised by a colleague of mine who is in- 
terested in a complete description of the solutions of the quadratic operator 
equation 

p2-2pp+(p2-g2)j„ = 0, p,gGM, with/-g2 >0 (1.2) 

where, in the following, denotes the unit n-dimensional matrix, and p is 
a Hermitian operator. When p is a density matrix, the equation ()1.2|) is the 
simplest generalization of the pure state condition = p. We note that 
Werner states ^1], and Horodecki states [0] satisfy such an equation. 

By itself, an equation as ()1.2() has nothing special. According to Cayley- 
Hamilton theorem, any finite n-dimensional matrix satisfies an n-degree 
polynomial equation, the characteristic polynomial. As we will see in the 
following, equation (|1.2)) describes a Hermitian operator whose spectrum 
has a maximum degeneracy. If the multiplicities of eigenvalues are denoted 
by k and n — k, respectively, the spectral decomposition of p is 

p = XlP + X2iIn-P) (1.3) 

where P is the projection onto the /c-dimensional subspace. Thus, a de- 
scription of the projection P on C" is equivalent to a description of the 
Grassmann manifold Gr{k,n). 
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Formally, the states of an n-level quantum system are described by den- 
sity matrices p that are positive, Hermitian, nuclear operators whose trace 
is normalized to unity 



p = {p > 0, p = p*, Tr p=l] 



(1.4) 



where * denotes adjoint, i.e. the complex conjugated transpose. 

If we denote hy D = (Ai, . . . A„) the diagonal matrix of p eigenvalues, 
elementary facts from the spectral theory of self-adjoint operators tell us that 
there exists a unitary operator U £ U{n), where U{n) € EndC^, generated 
by the eigenvectors of p such that 



with Y2i ^« = 1- III ths generic case, a Hermitian n x n matrix H is param- 
eterized by real parameters, number which coincides with the number 
of parameters entering the parameterization of an arbitrary unitary matrix 
U E U{n). Since the eigenvalues Aj are also independent, we infer from (|1.5)) 
that there are some constraints upon the form of U, i.e. the number of free 
parameters entering U is less than n^. It is well known that the constraints 
coming from the Hermitian character of an operator H may be translated 
to the request that its eigenvectors are defined up to an overall arbitrary 
phase. We choose the phases such that the first entry of each eigenvector is 
a non-negative number. This means that the first row entries of U are non- 
negative, i.e. U is parameterized by n(n — 1) real parameters, and in this 
way we conclude that a generic Hermitian matrix is parameterized by real 
parameters, as it should be. In fact this approach will be used everywhere 
in the paper: we start with a parameterization of U{n) and restrict to the 
appropriate subset of coordinates to describe the corresponding manifold. 

We made this digression because, especially in the physical literature, 
U, entering equation (|1.5|) . is considered an element of SU{n — 1), that 
evidently leads to an over parameterization. Thus when we have to do some 
symbolical, or even numerical calculations, as in or ^2j, we have to be 
more careful. We shall see later that such a J7 is a matrix realization of the 
flag manifold 



where [/(l)"" denotes the torus subgroup of U{n). By obtaining a full and 
explicit parameterization of X, we obtain, via formula ()1.5() . a parameteriza- 
tion of all finite-dimensional Hermitian operators whose spectra are simple. 



p = UDU* 



(1.5) 




(1.6) 
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Such a construction addresses an old fundamental question in the theory 
of measurement |15|.|S]. namely if it is possible to measure experimentally 
the 'variables' corresponding to an arbitrary Hermitian operator. Thus the 
first question to be solved is the finding of the 'variables' entering an Her- 
mitian operator. After that, the answer is simple: there does exist an ex- 
perimental embodiment for every Hermitian operator in finite-dimensional 
Hilbert space; see ^_UI for details concerning its realization. 

If /O is a density matrix, then p in equation ()1.5|) is parameterized by — 1 
parameters because of the trace condition. At the other extreme, there is 
the case when H is an one-dimensional projection, and then U entering (jl.Sj) 
is a matrix realization of the coset 



where G(l,n) denotes the simplest Grassmannian. Between these two ex- 
treme cases lie all the other spectral types. For example, the solutions of 
equation H1.2|) describe the most degenerate spectrum, which means that p 
has only two distinct non-zero eigenvalues. By description, we will under- 
stand the parameterization of the set of all unitary diagonalizing operators 
U in terms of a subset of parameters entering the parameterization of the 
group U{n), for a generic situation. In other words we are looking for a 
matrix realization of the relation (|1.3|) . i.e. of the projection P, and this 
is the place where the Grassmann manifolds enter the play; and we expect 
that in this case, stronger conditions on the form of U will be in force. In 
fact, in this paper, we provide a unified method for treating all the spectral 
types of Hermitian operators, by properly taking into account their spectral 
multiplicities. 

This problem is closely related to the description of isometrics between 
the Hilbert spaces C*' and C", 1 < k < n. The isometries are operators 
generated by n x A; or A; x n matrices whose columns, and respectively rows, 
are orthogonal, and in the following we show that there is a close relationship 
between these isometries and different matrix realizations of the coset spaces 
generated as in (1.1). The necessity of working with matrices that have 
orthogonal column vectors, or row vectors, became evident in the last years; 
see e.g. 0, H- 

The mathematical background necessary for obtaining such results are 
elementary facts from the spectral theory of self-adjoint operators and the 
theory of contraction operators, and a trivial lemma that we state here for 
the case of the n-dimensional unitary group G = U{n); for the general case 
see jSj. 
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Lemma 1. The coset relation 



X ^ U{n)/K 



(1.7) 



where K C U{n) is a subgroup, can he written as a matrix relation in the 
following form 



where Mn € U{n) is an arbitrary n x n unitary matrix, An G U{n) is a 
unitary matrix parameterized by a point of the coset X, and Bk ^ K is 
an arbitrary element of the subgroup K viewed as an element embedded in 



By using it we obtain new parameterizations of flag, Stiefel and Grass- 
mann manifolds. The paper is a sequel of our results concerning the factor- 
ization of unitary matrices j2j in terms of n complex vectors Vi € 5^*^^, i = 
1, . . . , n, where is the /c-dimensional sphere in C". We denote by M{n, k) 
the set of all n x A; complex matrices over C", and the main mathematical 
result of the paper is: 

Main Theorem. Let (t G M(n, k) be an n x k complex matrix that 
generates an isometry, i.e. C : C'^ — > C", C = Ik- Then, the matrix 
representation of the coset generated by the point (t is realized in terms of a 
unitary matrix A{n,k) that diagonalizes the projection Dt* , where Dt* is 
the defect operator associated to the isometry <t, under the form 



where Ci, i = 1, . . . ,k, are the (orthogonal) column vectors of <t. 

All the other results discussed in this work are a consequence of the 
above theorem. 

The organization of paper is as follows. In section 2 we show the close 
relationship between the isometries <t G M{n, k) and matrix realizations 
of the coset spaces. In section 3 we reformulate our results j2| in a more 
convenient form for the present applications, and we find the generic param- 
eterization of n-dimensional Hermitian operators. In section 4 we obtain our 
matrix parameterizations of Stiefel and Grassmann manifolds. The paper 
ends with Conclusion. 



Mn = An Bk 




U{n). 




(1.9) 
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2 Isometries 



In this section we show how the isometries C G M{n, k) can be used for 
the parameterization of various interesting manifolds. The main idea is 
that the isometries generated by k rows, or columns of an arbitrary n x n 
matrix, 1 < k < n, allow to define projection operators whose spectral 
decomposition provides the necessary tool in finding matrix representations 
of various interesting manifolds. For what follows, we need a few elementary 
notions from contraction operator theory; we use this theory since it has a 
powerful functional calculus that help us in doing explicit calculations. 

An operator T applying the Hilbert space !K in the Hilbert space !K' is 
a contraction if for any v G 'K, \\T v\\:yc' ^ l|f^l|:K) i-e. \\T\\ < 1, where ||T|| 
denotes the norm of T |13j . For any contraction we have T* T < I-ji and 
XT* < /jf/, where T* denotes the adjoint, that is defined by the relation 
(Tv,v') = {v,T*v'), V G "K^v' £ "K', and (•,•) is the usual inner product 
in !K, or, respectively, To any contraction T one associates two defect 
operators by the relations 

Dt = {I^ - T*T)^/\ Dt* = {I^' -TT*)^/^ (2.1) 

that are Hermitian operators in 'K and 'K' , respectively. They have the 
property 

TDt = Dt* T, T* Dt* = Dt T* (2.2) 

Finite-dimensional contractions can be seen as being generated by n x /c, or 
k X n matrices, and for definiteness we consider the first case, i.e. T has n 
rows and k columns and we denote it by C In the following we are interested 
in contractions of a special form, namely the isometries between C'^ and C". 
They are operators T : C'^ ^ C" that satisfy T*T = 4, k = l,...,n, i.e. 
the columns of T are orthogonal columns, and respectively, T : C" 
that satisfy TT* = J^,, when the k rows are orthogonal. With our choice, 
i.e. T is a n X /c matrix, we have the identification "K = and "K' = C". 

For an isometry £ generated by k columns, the relations (|2.H) give Dt = 
Dif = 0. However, in the following we preserve the notation for the defect 
operator Dt, to be implied that it is generated by a definite contraction, 
and from the first relation 1)2.2(1 we deduce that 

Dt* T = T Dt = = Dt* <t = Mt (2.3) 

i.e. all columns of <t are in the kernel of Dt* ■ In other words, the k column 
vectors are the eigenvectors that correspond to the eigenvalue A = 0, and 
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these eigenvectors are orthogonal. When Dt = 0, the other defect operator 
Dt* is a projection, i.e. a self-adjoint operator which satisfies Dt* = 
and we infer that rank Dt* = n — k; thus /„ — Dt* projects onto the 
eigenspace corresponding to the eigenvalue A = 1. Then, from the preceding 
relation, we get 

{In-DT*)<t = In(t=(t = X<L (2.4) 

i.e. the k column vectors of (t are the orthogonal eigenvectors of /„ — Dt* 
corresponding to the k eigenvalues A = 1, or in other words, rank (/„ — Dt* ) 
= k. 

For us the interesting objects are the unitary matrices that diagonalize 
the projections Dt* and In — Dt*- Let A{n, k) be the unitary matrix which 
diagonalizes the projection Dt*, then 

Ain,krDT*Ain,k)=(^ ^ ) (2.5) 

and from it we obtain a matrix representation of the projection Dt* under 
the form 

DT*=Ain,k)( / )Ain,kr (2.6) 



But we also have 
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Ain,knin-DT*)Ain,k)=( ^ ] (2.7) 



relation which can be written as 

k 

In - Dt* = A{n,k) ( ^ \ A{n,kr = ^ c, • c* = (2.8) 

^ ^ i=l 

where Ci, i = 1, . . . ,k, are the column vectors of the isometry <t. Equations 
(|2.6j) and (|2.8j) provide matrix representations for both the projections on the 
n — k-, and respectively /c-dimensional, subspaces of C". In our applications 
A{n, k) will be an n X n unitary matrix generated by a coset as in Lemma 
1, coset, which at its turn, is parametrized by the point C 

The above representation does not provide a full explicit form for the 
matrix A[n,k). We know only that its first k columns coincide with the <t 
columns. Thus an important problem is the completion of the matrix A{n, k) 
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with n — k columns without introducing new parameters, i.e. the next n — k 
columns must be determined by the first k columns. For doing that we 
need a parameterization of the first k columns and the most convenient one 
is to introduce generalized spherical coordinates. In the next section we 
show that the most important case is k = 1, the other cases being a direct 
consequence of it, which leads to a factorization of unitary matrices. In the 
same time, we reformulate our results from [2], concerning factorization of 
unitary matrices, that will provide the necessary tools for obtaining matrix 
realizations for symmetric manifolds. 

3 Factorization of unitary matrices 

The idea behind such a factorization comes from the following sequence 

rr, N Uin) U{n-l) U(2) 

U(n) ^ , V X —- f X • • • X —-4 X U(l) 

^ ' U{n-l) U{n-2) U{1) ^ ' 



^ ^^n-l X ^2n-3 x • • • X 5^ X (3.1) 

sequence that shows that each factor can be parameterized by an arbitrary 
point on the corresponding complex sphere. 

The factorization comes from the relation (3.1) which we write as 

g2n-i ^ u{n)/U{n-l) (3.2) 

Similarly to relation H1.8|) . we rewrite (3.2) as 

Mn = BlMl_^ (3.3) 

where £ U (n) is an arbitrary n x n matrix, £ U{n) is a special matrix 
generated by one of its column vectors, for example by its first column vector, 
which at its turn is parameterized by a point z € S"^""^, and 

= ( ; m!_. ) p-^' 

where M^-i G U{n — 1) is an arbitrary (n — 1) x (n — 1) matrix. For further 
details see [2j. By iteration of (3.4), we arrive at the form 

M„ = a0.i?i_i...i?ri (3.5) 



8 



where each B^_i^ is generated by a column vector, and here is the point where 
an exphcit parametrization enters. First we remind that a standard param- 
eterization of U{n) is given in terms of n(n — l)/2 angles, 9i G [0, 7r/2], i = 
1, . . . , n(n — l)/2, and n(n + 1)/2 phases, ipi G [0,2tt), i = 1, . . . , n(n + 1)/2, 
see e.g. [T], 

The matrix realization of formula (3.5), which is the main result in [2], 
is given by 

Theorem 1. Any element Mn € U{n) can be factored into an ordered prod- 
uct of n matrices of the following form 

Mn = Bl-Bi_^...Bl-^ (3.6) 

where 



uk _( h \ 

and Bn-k £ U{n — k) are special unitary matrices, each one generated by 
a single complex (n — k)- dimensional unit vector, bn-k & 52{n-fc)-i_ j^^^ 
example B^~^ = e*"^, where ip is an arbitrary phase. 
If Um £ 5*^™"^, m = 1, . . . , n, is parameterized by 

= (e^'^i cos 01, e*'^2 gin 0^ cos 02,..., e*'^™ sin 9i . . . sin Om-iY 
where t denotes transpose, the m columns of Bm are given by 



I e*^i cos Oi 

e*'^2 sin Qi cos 62 



\ e**^™ sin 61... sin 6*^-1 / 



and 



Vk+l 



dOk 



ViWi 



Ok-i = ^/2), 



1, m-1. 



where in the above formula one calculates first the derivative and afterwords 
the restriction to 7r/2. 
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The essential point of the algorithm for getting (|3.6|) was the use of 
contraction operator theory to construct all the eigenvectors of a self-adjoint 
projector, generated by the first column vector of B^, eigenvectors that are 
easy provided by the above theorem, if the first vector generating them is 
parameterized in spherical coordinates. 

Taking into account that by multiplication at left of an arbitrary matrix 
from U{n) by a diagonal phase matrix dn = (e*"^^, . . . , e*'^"), the first row is 
multiplied by e"^"^^ , the second by e^'^'^, etc., and the last one by e^'^" we can 
write = d^B^ where the first column of B^ £ SO{n) has non-negative 
entries. 

For our aims we need an explicit parameterization of Bn-k, k = 0, . . . , n— 
1, and we choose the n generating vectors as follows 

Vn = cos oi, 6*"^^ sin oi cos 02, . . . , , e*"" sin oi . . . sin a„_i)* 

Un-i = (e*'^^ cos 61, e*^2 sin bi cos 62, ... , e*^""i sin bi . . . sin 6„_2)* 
: (3.7) 

2/2 = (e*^^ cos zi,e^^^ sin zi)* 

yi = e*"^! 

The projection operator /„ — Dt* entering formula (|2.8|) is a self-adjoint 
operator, and for such an operator its eigenvectors q are defined up to 
an overall phase. We choose the phases such that the first entry of each 
eigenvector is a nonnegative number. With the generating vectors of the 
form (|H.7)) , the parameterization of the matrix , given by Theorem 1 , 
is such that its first row entries have the form 

mil = e*"^ cos ai, mi2 = — e*("^'''^i) cos bi sin ai, 

mi3 = e*("i+^i+'>'i) cos ci sin oi sin hi, (3.8) 

min = (-1)""^ gi(ai+-+i^i) sin ai . . . sin zi 

and if we want that these matrix elements should be nonnegative we have 
to take 

«! = 0, /3i = 7i = ■ ■ • = wi = vr (3.9) 

In the following we shall use in relations 1)3. 7|) these constraints, and so 
we remove the first phase of each of the vectors yk and change the numbering 
as ai — > ai-i, i = 2, . . . , n, /3i — > Pi-i, i = 2, . . . ,n — 1, etc., and now each 
vector yk is parameterized by 2{k — 1) parameters, k = 2, . . . ,n, i.e. an 
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equal number of phases and angles. In this way the last vector disappears, 
such that the relation H3.6|) has the following form 



Tln = B^^-Bl_,...Bl^-^ (3.10) 
The last matrix is the matrix realization of the flag manifold 

. U(n) U(n) 

Fl{n) ^ — ^ —H- 3.11 

^ ' U{1) X U{1) X ••• X U{1) U{lY ^ ' 

it depends on n{n — 1) parameters and it is the most general form of a 
unitary matrix that diagonalizes an n-dimensional Hermitian operator H 
whose all the eigenvalues are simple. In the same time Fl{n) is the 'natural' 
manifold on which one can define a symplectic structure, and an interesting 
problem would be to find its explicit form. 

Putting together the information contained in Theorem 1 and formulas 
(|.S.in|) and (|3.11|) we have the following 

Corollary 1. The matrix 971^ describes the spectral decomposition of a 
finite- dimensional Hermitian operator H whose eigenvalues are simple, i.e. 
they satisfy a relation as 

Ai > As > • • • > An 

and the operator H is written in the form 

n 

H = Y,\uiul (3.12) 

where Ui, i = l,...,n, are the column vectors of the matrix \3.1CI\) . If 
H > and Tr H = /i E M^, then Aj, entering formula (3.12), could he 
parameterized as 

Xi = h cos 6*1, X2 = h sin 9i cos 6*2, ... , Xn = h sin 6*1 .. . sin 0„_(3.13) 

where 6i € [0, 7r/2], i = 1, . . . , n — 1, are arbitrary angles. If H is a density 
operator, h = 1 in the above relation. 

If H is not positive definite, let suppose that its first p eigenvalues are 
positive, and n — p are negative. If Tr H = 0, and h = TrH\.yQ = 
—TrHx-^Q G M*^, then a parameterization of eigenvalues is given by 

O O O / \ 

Xi = h cos 9i, X2 = h sin 9i cos 02, ... , Xp = h sin 9i . . . sin 0p_i(3.14) 
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Ap+1 = —h cos^ 9p, . . . , Xn = —h sin^ 6p . . . sin^ On-i (3.15) 
// Tr H = h ^ the parameterization of eigenvalues is given by 
Ai = \h\ cosh^ 6, A2 = \h\ cosh^ 6 cos^ 6*1, . . . , 

Xp = \h\ cosh^ 9 sin^ 6*1 .. . sin^ 6lp_2 (3.16) 

Ap+i = —\h\ sinh^ 0, . . . , A„ = — sinh^ sin^ ^p-i • • • sin^ (3-17) 

if h > 0, where 9 E M^, 9i S [0, 7r/2], z = 1, . . . , n — 2, are arbitrary angles, 
and by a similar formula in which one interchanges ch^9 ^ sh'^9, ifh<0. 

In this way, Corollary 1 gives a simple and explicit parameterization 
of all generic finite-dimensional Hermitian operators. 

4 Matrix realizations of Stiefel and Grassmann 
manifolds 

Looking at relation 1)3. 11(1 we consider that the next simpler manifold is the 
Stiefel manifold that we define as the coset space 

U{n) _ Fl{n) 

- U{lfxU{n-k) - Fl{n-k) ^ 

instead of the usual definition 

Both forms are similar, only the number of parameters entering them is 
different. According to Lemma 1 we can write the first relation (|4.1() as a 
matrix relation 

M„ = ^(„, (4.3) 

where M„ € U{n), Mn-k £ U{n — k) and A{n, k) € U{n) is the matrix real- 
ization of the Stiefel manifold St{k^ n) parameterized by a point represented 
by a n X complex matrix, with the first row entries non-negative numbers. 
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Alternatively, we may consider that M„ G Fl{n), Mn-k Fl{n — k). Look- 
ing at relation ()2.8|) we see that A{n, k) is the same object in both relations 
H2.8|) and (4.3). By consequence we deduce from 1)3. 10() that 

A{n,k) = BlBi_,...B^-l^^ (4.4) 

is the matrix realization of the Stiefel manifold ()4.1() . Thus the following 
holds 

Theorem 2. The matrix representation of the Stiefel manifold St{k,n) as 
defined by (|4.1j) is obtained from the parameterization (|3.10|) by taking zero 
all the parameters entering U{n — k), i.e. U{n — k) = In-k such that 

A{n,k) = Bl.Bi_,...B^^-_l+, (4.5) 

By using elementary facts from spectral theory of self-adjoint operators we 
find that the matrix representation of the projection In — Dt* writes as 

In - Dt* = St{k, n) =cicl + --- + Ckcl (4.6) 

where Ci, i = 1, . . . ,k are the first k column vectors of (|4.5() . As long as the 
matrix (|4.5|) is parameterized by di = — k — {n — k)"^ = k{2n — k — 1) real 
parameters, we can choose any k column vectors of (j4.5j) in formula (4.6) 
and we can make this choice in (^) modes. If the form (|4.2() is used, then 

A{n,k)=B^,.Bl_,...B'^ZU, (4.7) 

where Bf^^p, p = 0, . . . ,k — l are the first k matrices entering equation (|3.6|) . 
i.e. the complex vectors generating (4.7) are such as in the written form 
(|3.7|) . comprising n(n + l)/2 phases. The difference between the forms (4.5) 
and (4.7) consists in the number of real parameters entering them: (4.5) is 
parameterized by di parameters, while (4.7) depends on d2 = n'^ — {n — k)'^ = 
k{2n — k) parameters. 

Remark. Similarly to Fl{n), on the Stiefel manifold (|4.1j) one can define 
a symplectic structure. 

As a consequence of the above theorem, we have the following 

Corollary 2. The matrix (|4.5|) describes the spectral decomposition of a 
finite- dimensional Hermitian operator H that has a degenerate eigenvalue 
of multiplicity k, Xi = ■ ■ ■ = Xk, and all the other eigenvalues are simple. In 
this case H writes as 

k n—k 

i7 = Ai ^ Uj n* + ^ Afc+i Uj+fc u*^^ (4.8) 
i=i i=i 
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where Ui, i = 1, . . . ,n, are the column vectors of the matrix H4.5|l . If H is a 
density operator, then 

r\ 

Ai = cos 01, Xk+i = sin 0i cos 02, ■ ■ ■ , Xn = sin 0i sin 02 ■ ■ ■ sin 6'n-fc-i 

where 0i € [0, vr/2], i = 1, . . . ,n — k — 1 are arbitrary angles. If H is not 
positive definite, the eigenvalues are parameterized by similar formulas, as 
in Corollary 1. 

The Grassmann manifold, Gr{k, n), is defined as the set of all fc-dimension- 
onal subspaces of C", and a main problem is to have a simple description 
of this variety. From what was said before, evidently there is a one-to-one 
correspondence between A:-dimensional subspaces, fc-dimensional projections 
and a unitary operator like A{n,k). The first problem is to determine the 
number of independent parameters entering A(n, k) for this case. Taking 
into account that k- and n — fc-dimensional subspaces are transformed into 
k- and respectively n — A:-dimensional subspaces under the action of matrices 
from U{k) C U{n) and respectively U{n — k) C U{n), Gr(k,n) is viewed as 
the coset space 



Gr(kn)^ ^ - ^ (4 9) 

^''~Uik)xUin-k) Fl{k) X Fl{n - k) ^ ^ 

and by consequence the real dimension of the Grassmann manifold Gr(k, n) 
equals 



d = n^ -k^ - {n-kf = 2 k{n - k) (4.10) 

Like the case of Stiefel manifolds, we want first to obtain a parame- 
terization of A[n, k) which is equivalent with finding a parameterization of 
Grassmannians. With that end in view we rewrite relation (4.9) into the 
form 

where M„ S U{n) is an arbitrary matrix from U{n) and € U{k) and 
respectively Gn-k € U{n — k). Such a relation is always possible since any 
group operation is uniquely written as the product of an element in the 
subgroup with an element of the coset. From that relation we infer that 
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the number of independent parameters entering A{n, k) is the same as the 
dimension d of the Grassmann manifold Gr{k,n). 

If we look at the relations H3.1U() or H4.5() and consider one of them for 
the case A; = 1 we observe that B!^ is equal to A{n, 1). Thus we obtain a 
matrix representation for 

In-P = Gr(l,n) = Ain, 1) J ^ ) A*{n, 1) (4.12) 

i.e. the simplest Grassmannian, a result already known. The preceding 
relation can be written in an equivalent form as 

Gr{l,n) = vi-vl (4.13) 

where vi is the vector that generates B^, i.e. the vector y„ from relation 
(|3.7j) . Thus A{n,l) can be obtained from (|3.1()|) by taking all the phases 
and angles entering U{n — 1) equal to zero, i.e. U{n — 1) = In-i, and here 
we give its explicit form 



/ 1 



A{n,l) 



(4.14) 



/ cos a\ — sin a\ . . . \ 

sin a\ cos a2 cos oi cos 02 — sin 02 . . . 



- sin a„_i 

\ sin ai . . . sin a„_i ... ... cos a„_i / 

where we factored out the diagonal phase matrix. This formula provides an 
explicit form of the vectors Vi that are generated by Theorem 1. 

In the following we show that A{n, k) can be obtained in a similar way. 
Taking into account the form of the projection operator In — P and the 
dimensions di and d for Stiefel and Grassmann manifolds respectively, we 
infer that the Grassmann manifold is a special case of a Stiefel manifold. 
Our problem now is to find those constraints which lead to the correct 
parameterization of A{n, k) for Grassmannians. 

In order to view which are the constraints, let us consider the case k = 2. 
By taking into account that the first column of A{n, 2) coincides with the 
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first column of the matrix ()4.14|) , we infer from H4.10() that the parameteriza- 
tion of the second colmim is given in terms of 2(n — 3) new real parameters. 
This means that we have to remove an angle and a phase from the vector 
Un-i, equation ()3.7|) . A convenient choice is to take equal to zero the last 
angle and phase, i.e. 6„_2 = Pn-2 = 0. This choice induces the following 
form for the matrix B^_i 

(I \ 

SLl(ftn-2 = l3n-2 = 0) = B];\ = Bn-2 (4.15) 

V I I 

where Bn-2 is generated as in Theorem 1 by the vector 

y'n~2 = (cos 61, e*^^ sin 61 cos 62, ... , 6*^"'=* sin bi . . . sin 6„_3)* (4.16) 

It is easily seen that this structure preserves from /c — > /c + 1 such that 

(Ik \ 
b':\= B„_2fc , fc=l,...,[^] (4.17) 
V 4 / 

where [a] denotes the integer part of a, and -Bn-2fc; A; = 1, . . . , [|], are gen- 
erated by the vectors 

wi = (cos oi, e*"i sin ai cos 02, . . . , , 6*°"-^ sin ai . . . sin an-if 

W2 = (cos 61, e"^^^ sin 61 cos 62, • • • , e*^"~^ sin 61 . . . sin hn-zY 

Wp = (cos /i, e"^i sin li cos h,-- - , e*'^"-(2p-i) sin /i . . . sin ^„_(2p-i))* (4.18) 

^[n/2] = (cos zije*"^! sinzi cos 22, e*'^^ sin 2:2)*; for n odd 

i(;[„/2] = (cos zi^e^^^ sin zi)*, for n even 

With the above notation the following holds 

Theorem 3. The unitary matrix A{n, k) € U{n) entering the matrix rep- 
resentation of Grassmann manifold has the form 

A{n, k) = B^^- Bl;\ . . . ^ = 1^ • • • J^] (4-19) 

and the matrix representation of the projection onto the k-dimensional sub- 
space of C" is 

In-DT* =Gr{k,n)=A{n,k)(^^^ °^^(n,A;)* (4.20) 
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An equivalent description is given by 

k 

Gr{k,n) =^Ui-u* (4.21) 

1 

where Ui, i = 1, . . . ,k, are the first k column vectors of matrix (4.19 ). 

Corollary 3. If H is a Hermitian operator that has only two eigenvalues 
whose multiplicities are k and n — k then 

k n 

H = XiYl + ^2 5^ u,u* (4.22) 

i=l i=k+l 

where Ui, i = 1, . . . ,n, are the column vectors of the matrix A{n, k), equation 
(4.19). If H is a density matrix operator, then Ai = cos^ ^ and A2 = sin^O, 
where 9 S [0,7r/2] is an arbitrary angle. With this choice, formula (4.22) 
describes all the density matrices that are solutions of the equation (1.2). 

When H is not positive definite, and Tr H = h, the parameterization 
of eigenvalues is: Xi = h cosh'^ and A2 = —h sinh^ 0, Ai > IA2I, and 
Ai = \h\ sinh^ 9 and A2 = —\h\ cosh^ 9, if h < 0. . 

In this way we succeeded to find an explicit parameterization for the 
finite-dimensional Hermitian operators, in three cases: the generic case, 
when all the eigenvalues are different, an intermediate case when the first 
k eigenvalues are degenerated, the other eigenvalues being simple, and the 
most "degenerated" case, when there exist only two (different) eigenvalues. 
Matrix realization of other flag manifolds, as e.g. 

Fl(v) 

Flikuk2,...,ki,;n)^ > — — (4.23) 

Fl{ki) X Fl(k2) X • • • X Fl{ki) 

with Yli f^i ^ that describe other spectral multiplicities, can be found in 
a similar way. Assembling together all the previous information we have 

Corollary 3. The flag matrix ()3.1U() encodes all the possible spectral de- 
compositions of finite- dimensional Hermitian operators, acting on C", i.e. 
any other particular unitary matrix A{n, k) that describes a given spectral 
multiplicity of an Hermitian operator can be obtained from it by properly 
restricting the number of parameters entering (|3.10() . 

Example. We consider that a simple example of a matrix realization of 
a Grassmann manifold will help, and we consider the case n = 8 and k = 4; 
then A{8, 4) has the form: 
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^(8,4) = 



cos ai 



\ e*"'' sinai . . . sin 07 



— sin ai 



sinai 














Sin 07 



e*"'' cos 07 / 



/ 1 

— cos 61 

e'^i sin bi cos 62 



V 





sin bi 
e*'^! cos bi cos 62 






— e*'^! sin 60 



e*'^^ sin 5^ . . . sin 65 










cos 65 

1 



/ 12 











— COS Ci 

o«7l sin d COS C2 





sin ci 

COS Ci COS C2 









sm C2 



e*'''^ sin ci sin C2 cos C3 6'"^^ cos ci sin C2 cos C3 e''^^ cos C2 cos C3 



=*72 



,«72 



e''''^ sin ci sin C2 sin C3 




e*'''^ cos ci sin C2 sin C3 




e''''^ sinc2COSC3 








-e*'^^ sinc3 
e''>'3 cos C3 











I2 



1 ^3 








\ 





— cos di 


sin di 








e'^^ sin di 


6**^1 cos di 





V 








^3 / 



(4.24) 



and the matrix representation of Gr{k, 8) is 

i=k 

Gr{k,8) = ^UiU*, for A; = 1,...,4 (4.25) 

i=l 

where Ui, i = 1, . . . ,k are the first k columns of matrix (4.24). In the same 
time formula (4.24) provides us all the Grassmannians Gr{k,n) for n < 7. 
For example A{k, 7) is obtained from (4.24) by putting 07 = 65 = C3 = di = 
a7 = /35 = 73 = (5i = and by deleting the last row and column, and so on. 
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5 Conclusion 



In this paper we have obtained a constructive parameterization of all finite- 
level Hermitian operators. Further, this construction is recursive: if we 
have a parameterization of Tin, is obtained by embedding Tin into 

9Jln+2i and by multiplication at left by an appropriate matrix of the form 
H4.14() . with n — > n + 2. We have shown that the unitary matrices which 
diagonalize the Hermitian operators are subsets of the flag unitary matrix 
Tin whose explicit form is given by formula (|3.1U|) . When the eigenvalues 
are simple, the generic form of the unitary operators is Tin', when there 
is a fe-fold degeneracy, the corresponding unitary matrices are in the set 
St{k,n). If there is a maximum degeneracy, i.e. only two eigenvalues, with 
multiplicities k and, respectively n — k, are distinct, the unitaries are in 
the Grassmannian G{k,n). Our explicit construction was done only for the 
n-dimensional unitary group U{n), but it is evident that the same approach 
works in the case of any compact group. For example, taking zero all the 
phases entering U{n) and doing similar calculations, one gets results for 
SO{n), and so on. 

We consider that the above parameterization will be useful for doing 
calculation, especially for problems where one has to make an optimization 
over a set of parameters, e.g. for the characterization of entanglement. 

Taking also into account our previous results j^, which state that any 
unitary matrix entering Tin is an ordered product of n — 1 diagonal phase 
matrices and n(n — l)/2 two-dimensional rotations, one can make use of 
the device designed by Reck et al ^01 to give an operational meaning in 
the real world to any finite-level Hermitian operator. This means that all 
the Hermitian matrices could be experimentally implemented and, by conse- 
quence, could be measured. Such multi-states devices will find applications 
in quantum information processing and quantum computation. 

In the same time, we obtained a new and simple analytic representa- 
tion of Stiefel and Grassmann manifolds, representation that is essentially 
contained in a unitary matrix A{n, k), which can be easy stored into a com- 
puter, and problems similar to those encountered in |3, or ^| will be easier 
to solve. 

Our results show also that there are two new symmetric spaces Tin, 
and St{k,n), defined as in (|4.1j) . on which one can implement a symplectic 
structure, and an interesting problem would be the finding of the relationship 
between a generic Hermitian operator on these symmetric spaces and the 
Laplace-Beltrami operator on the corresponding manifold. Some of these 
problems will be considered elsewhere. 
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